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(oN KO DUG lL LON 


The phenomena of physical Sputtering have p2en 
extensively studied Since it was first mentioned in 
Philosophical Transactions by Grove in 1852. Sputtering is 
the ejection of atomic material from a target resulting fron 
bombardment by an energetic particle. Most scientific 
investigation has been restricted to ionic bombardment 
because the ion can easily be accelerated to any desired 
velocity. Sputtering has defied complete physical 
description because of complicated mechanisms, difficulty in 
treating the behavior of a many-body system, and an 
incomplete understanding of the interacting forces and 


potentials in crystal structures. 


Sarly experiments typically determined the "sputtering 
eave 0 ' of Single and polycrystalline materials. The 
Paper recraing ratio” (interchangeable with sputtering yield - 
Y) is defined as the number of target atoms ejected from the 


surface per incident ion. 


Pieeeio2copeeKingdon and Langmuir { 1] bombarded thoriatced 


j4 


tungsten with various ions ina glow discharge tube. alee 
waS a special case of sputtering since the thin surface film 
of thorium on a tungsten substrate was sputtered rather than 
the eungstenm. itself. The results of this experiment 


qualitatively suggested an ejection mechanism. 


A few years later, Blechschmidt and von Hipple [2,3] 
proposed a theory which described sputtering as an 
evaporation of surface atoms (an emission mechanism). 


Earlier, von Hipple [4] had found by spectroscopic analysis 
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that at least some sputtered atoms were in an excited state. 
The sputtering theory showed that atoms in the region of 
mupact could rapidly acquire thermal energy if the kinetic 
energy of the incident ion was converted to thermal energy 
at the target surface. If an atom acquired a sufficient 
amount of thermal energy, it would then evaporate from the 


surface, some atoms evaporating while in an excitad state. 


Direct application of statistical methods to sputtering 
was made by Harrison {5] who envisioned the interaction of 
mom dilStribution Eunctions (the crystal lattice and the ion 
beam). These nodels, based on statistical methods, 


implicitly accept ejection-type mechanisms. 


One of the most important contributions to the study of 
Sputtering was made by Wehner [6] in 1953. He -s discounted 
the evaporation theory and presented strong evidence for a 
Momentum transfer process. Shortly after Wehner's findings 
were reported, Henschke [7] proposed a theory of sputtering 
based solely on classical theory, treating normal and 
oblique incidence sputtering separately. Hash CONG 2p not 
normal incidence sputtering required many-body collisions in 
which the ion was eventually reflected outward to sputter 
surface atoms. This theory was plausible for oblique 
incidence sputtering; however, the case for normal incidence 
Sputtering required that the ion be r2rlected inside the 


Sevsctal. 


Wehner and Rosenberg [8] reported that in the sputtering 
Of polycrystalline metals at low ion e@nergy (normal 
incidence) more material is ejected obliquely than in a 
direction normal to the target surface. This is due to the 
fact that less than a full (180°) overall direction change 
of momentum is more favorable for sputtering because this 
requires fewer number of successive collisions for aton 


Sqyeccion. In 1960, Wehner and Anderson {9]J reported 
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Sputtering patterns from single crystal targets which showed 
an anisotropic spread of energy from a collision center and 
atoms which are preferentially ejected in close-packed 


Teck ions. 


When an atom collides with another which has a nigher 
mass, the scattering angle can vary between 0° and 180°. [In 
a collision with a lower mass atom the scattering is 
confined to between O09 and 90° (thus backward reflection 
Cannot occur). In low-energy sputtering of multi-~eiement 
targets, Winters and Sigmund [10] showed theoretically that 
festgmiticant contribution to sputtering, called “single 
reflective collisions", arises from the reflection of lower 
mass atoms from underlying heavier atoms. For a detailed 


discussion of momentum reversal, see Harrison and Magnuson 


(tel 


mime ooe wewehnere ect al [ 12] published a compilation of 
beret umgeyreld> for metals and semiconductors in the 
100-600 eV range. In 1966, Wehner, Anderson, and KenKnight 
{ 13] published the results of a study that covered almost a 
decade of research in low-energy sputtering. A few of the 
numerous topics discussed included more accurate and 
improved techniques for analyzing the sputtered mass, the 
temperature dependence of sputtering yields of varioious 
metals bombarded by atomic and polyatomic ions in the 2.8 to 
10 keV range, the effect of oxygen on ion etection patterns, 
and the spatial distribution of ejected atoms in metals and 


semiconductorcs. 


A few low-energy sputtering experiments [{14, 15, 16, 17, 
18], which lend themselves to direct comparison with these 
Simulations, investigate the behavior of the sputtering 
Piao LUunCe#on Of 10h energy, lattice orientation, 


target material, and ion angle of incidence. 
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Paralleling the experimental investigation Of tne 
Sputtering yield and the development of mass spectrographic 
analysis was the attempt to model the _ genevering mechanism 
through the use of computer Simulation. In 1960, Gibson, 
Goland, Milgram, and Vineyard {19], the first investigators 
to use computer Simulation as a means of studying radiation 
damage, set forth the criteria that must be satisfied before 
the Simulation method could be applied to. crystals. 
Reference {19] gives a very detailed description of the 
Simulation method. Basically, they built a computer nodel 
to represent metallic copper and studied radiation damage 
events at low and moderate ranges of energ7 up to 400 eV. 
In their model, one atom in the stationary lattice? was given 
Speco tbary Kinetic energy and direction of motion, as 
though it had been struck by a bombarding particle, and was 
allowed to interact with the remaining atoms in th2 Lattices. 
This resulted in cascades of independent binary collisions. 
It is well established that sputtering of solids by ion 
bombardment is the result Oe successive independent 
approximately binary collisions somewhat as if these 
collisions take place with or between separate gas atoms. 
Low-energy ion scattering studies [20] have shown that the 
approximately binary scattering behavior still holds at 
bombarding energies as low as those approaching the 


sputtering thresholds. 


Mmiotseharrison, Levy, Johnson, and EEfron { 21] used a 
computer to simulate the bombardment of face-centered-cubic 
(fcc) copper single crystals by an argon ion. From their 
Studies, the sputtering mechanisms in the fcc crystal were 
isolated and identified. They reported, for ion energies 
Hess than 10 kev, the sputtering process occurred 
predominately within three atomic layers of the surface. 
The computer Simulation included only a repulsive force 


between the atoms. The equations of motion for the 
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many-body system were computed using an "average-force" 
algorithm rather than a "central-difference" algorithm, see 
Poreeson, Gay, and Efron [22]. In 1972, Harrison, M‘oore, 
and Holcombe [23] in continuing efforts to develop a more 
mB aiestaG Computer model incorporated a Cu-Cu potential 
function which included an attractive potential, surface 
tayersrelaxation, and Surface binding energy for tne (100), 
(110), and (111) copper planes. They reported improved 
agreement between experimental and simulation results, but 
the original interpretation of the low-energy sputtering 


mechanisms in fcc copper crystals remained unchanged. 


Computer Simulation of physical processes is still ina 
very infant state because the available processing speed of 
existing computers severely limits the sophistication of the 
models. However, past investigations were quite successful 
in identifying the mechanisns which occur when energetic 
particles interact with crystal atoms. Pie atl) wyl 2 o 
observe these mechanisms is an advantage peculiar to the 
computer Simulation because each crystal atom is identified 
by a number for the mathematical calculation, and its 
complete track can be labeled and recorded. The tracks of 
selected atoms can be examined to determine the dynamics of 


Peemaeoansuwnich Exit the front surface. 
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Wine oO haAkem GoubCcrll VES 


Peau LAR DLISPERSION OF ISOTOPES 


Recently, Wehner [24] reported that the isotope ratio in 
metal deposits obtained by sputtering a flat target at low 
10n energy (100eV) under normal incidence is a function of 
the ejection angle. He investigated the bombardment of a 
Cu(63-65) target with 100 eV Hg ions . He reported that the 
lighter isotopes are enriched in a direction normal to the 
target surface. This research was a computer simulation 
investigation of Wehner's studies of angular dispersion of 


1sotop2s in a monocrystalline copper Lattice. 


Peel lQuULS SPMASH EFFECT 


In preliminary NPS simulation studies, after analysis of 
data, a predominance of sputtered atoms occurred at certain 
distances from the center of the lattice representatives area 
Miwiniehiecme LON L@pacts Occur. Figures 1 and 2 [2521 depict 
the top layer of a copper fee crystal respectively bombarded 
by neon and gold ions at 2 keV. The numbers represent the 


number of times that the atom was sputtered in 365 events 


(Shots) . The dark triangle in the center represents the 
target impact area. All 36 events are directed into this 
representative area. The darkened elipses highlight che 


atoms most often sputtered and suggest a ring outward from 


the center (especially for the Cu-Au interaction). 
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Macklin and Metaxas [26] experimentally investigated the 
splashing of liguid drops on deep and shallow liquid layers. 
They published high speed photographs which showed liguid 
splashes that might be analogous to surface atoms be2ing 


ejected from the surface and a few eventually sputtering. 


This research was an investigation of "ring formation" 
or "liguid splash" effect of sputtered atoms. Determination 
was made of the probability that a congregate of atoms at a 


given distance from the impact area reference will sputter. 


Secu erGY CCNTAINMNENT 


AS with any computer simulation, a trade-off must be 
made between the accuracy and detail desired and the size 
and time requirements imposed by computer availability. 
Because of this constraint a limit exists on the size of the 
monocrystalline lattice used in the Simulation. This 
research was an investigation to determine the size lattice 


required to contain a low energy event. 


Pee LON VATOM MASS RATIO 


In July of 1976, Eer Nisse [27], reported experinental 
confirmation of a monotonic dependence of Y with the atomic 
number. He investigated the yield of a gold target impacted 
by 40 keV ions up to the atomic number of gold. This 
monotonic dependence is predicted by a sputtering theory of 
Sigmund [28]. However, computer Simulation data of previous 


NPS studies (29] indicate that the monotonic dependence of 
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the sputtering ratio on the ion mass ceases to be valid for 
lon/atom mass ratios greater than one. This research was an 
investigation of the sputtering ratio versus lon atomic 


number for a wide range of ion/atom ratios. 
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ieee! GALTON MODEL MODIFICATIONS 


The kEasic simulation model used in this research has 
been developed by the NPS group over the last decade [21, 
Pere 31, 32 ]- A brief description of this model and the 
alterations made to it for this research follow. POE 
detailed discussions of the simulation method, the read2r is 
@erected to references 19, 20, 21, and 23. A detailed 
description of the mathenatical model has been published 
n22 }. 


In the development of the computer model, the lattice 
unit (LU) proves to be a convenient unit of length, and is 
used extensively. One LU=a,/2, where a, is the cubic 
lattice constant, the atomic spacing as determined by x-ray 


S@meysceditoOogtaphic Studies. For copper, a, =3.615 angstroms. 


Certain properties of the program and results are best 
considered in the framework of a specific run. Each run is 
begun by constructing a target microcrystallite whose sites 
represent the equilibrium positions of copper atoms near the 
Slieieecomen a ECC Crystal. The target atom - target atom 
interaction is represented by the Born-Mayer Gibson type II 


(Gibson-II) repulsive potential function. 


Each ion of the beam is represented by a single neutral 
atom of the desired element [30]. The ion is placed at the 
DoEnemmOonme ltS trajectory such that its velocity vector 
determines the point of impact. The term bullet and ion are 
used interchangeably throughout this report for the incident 


particle to avoid confusing it with target atoms. 
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The bullet-target potential function was always either 
the Gibson-II or that determined by the Kinetic Secondary 
Electron (rs 2-5) Potential [30]. These two potential 
Functions seem to represent a reasonable range of "hard" and 


"soft" potential functions. 


As reported by Harrison, Moore, and Holcombe [23], the 
attractive forces between the atoms in the crystal are so 
small, when compared to the forces which occur in aven 
moderately energetic collisions, that they may be safely 
neglected in sputtering Simulations. Aliso, in order to use 


the smaller reduced target areas, the lions must be normally 


incident. Even though the computer program still retained 
the capability =O utilize attractive potentials and 
obliguely incident ions, they were not used in this 


investigation. 


A. MICROCRYSTALLITE SIZE 


Initial discussion to establish research objectives 
revealed the advantages gained from uSing a symmetrical 
Mitttoo that Could be fEolded and rotated so that atoms 
Sputtered by ions incident at all possible lmpact points 
would beccme available for analysis. A lattice with each 


layer a square and naving an odd number of atoms along each 


boundary was needed. The lattice generation program as 
previously used CS) pili handle aS change with no 
alterations. 


Figures 3, 4, and 5 show a roughly isometric projection 
of the labeled atoms in a 9x4x9 (100) orientation, a 9x5x9 
(110) orientation, and a 9x4x9 (111) orientation of a fcc 
Microcrystallite respectively. Aicialisy “ofl (cue Me js)elora| (atoye ation eur 


first used by Harrison and Delaplain { 31], with the distance 
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between layers in the Y direction exaggerated, and the atoms 
represented by 20° ellipses has proved to be a valuable tool 
for analyzing the computer programs and for displaying data. 
The first number in each lattice is atom number two, with 


atom number one being the bullet. 


From previous studies, it was clear that at the low 
energy levels considered the sputtering mechanisms were 
confined to the surface layers [21]. For this reason, the 
(100) and (1117) orientation microcrystallites were only four 
layers deep. The (110) orientation lattice generation 
program constructs a lattice which has Significantly fewer 
atoms and its layers are closer together. In an attempt to 
increase the number of atoms in the lattice while keeping 
the required computer tims below that of the (100) and (111) 
orientation programs, the (110) orientation microcrystallite 


was set at six layers deep. 


Dee cRLOLCAL BOUNDARIES 


The previous Simulation program contained a 
microcrystallite with established boundaries that were not 
consistent with the development of the potential functions. 
Slnce the potential Luncte Lon is truncated at the 
nearest-neighbor separation (ROZ), all boundaries (except 
the top surface) were set at ROE. This is justified by the 
fact that an atom outside these boundaries is farther from 
bnew cEYstal than the potential truncation radius. The upper 
surface boundary was set by starting a small distance from 
the ie ealger= surface, and increasing the distance in 
succeeding runs until the total yield was independent of “he 
distance. For computational efficiency this distance should 


be set as small as possible. 
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Peaqacre 6G showWwS a drawing of the microcrystallite in 
solid lines and the plane boundaries in dashed lines. The 
boundaries which are defined in the Glossary (Appendix C) 
are shown as planes at the stance E@ tne 


microcrystallite surface cited above. 


Se ESOTOPE GENERATION 


In order to study the angular dispersion of sputtered 
lsotopes, a method for generating a microcrystallits with 
the correct isotope ratio and masses was needed. The method 


used assigned an additional variable to each atom that 


identified the atom "type", The atoms were classified 
according tc two types; type 1- least abundant, heavier 
isotope and type 2 - most abundant, lighter isotope. 


Appendix A is a detailed discussion of how the correct 
isotope ratio was obtained and each atom properly identified 


With the correct mass. 


D. TARGET IMPACT AREAS 


A representative area must contain all possible impact 
Pemmts Upon a particular surface, that is any point in the 
Surface must be projectable Lie a point of the 
representative area. Figure 7 shows the intrinsic geometry 
of the three fcc orientations and the complete sets of 
impact points for each representative area. Figure 8 shows 
the reduced sets of impact points from which th2 complete 
sets are generated by reflection and rotation. Figures 9, 
10, and 11 show an expanded view of the target impact areas 
and the actual impact points used. The individual impact 


points are those used in previous Simulations for normal 


‘ 
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incidence. fiirs —ottplvtveatton, Which greatly reduces the 
computer running time, is only possible when the ion is 


normally incident. 


Eee GALTICE FOLDING 


Myonrdcr tosinvestigate the ring formation or "liquid 
splash" effect, it was necessary to develop a means of 
folding the results obtained from the reduced target area. 
This requirement was satisfied by writing three subroutines; 
MmEoOEoA, nrOLDB, and KFOLDC. Appendix B gives a detailed 
@escription of the "folding" logic and the geometry of each 
Orientation. Figure 12 shows how the (100) and = (110) 
Orientation reduced representative areas were folded to 
obtain the complete representative areas. Figure 13 shows 
how the (111) orientation reduced representative area was 
folded and summed to obtain the complete representative 


area. 


fen cOn SeUTTERING PROBABILITY 


A means of taking the folded data and calculating a 
conditional probabilty that an atom would sputter from a 
particular ring was needed. This was accomplished by 
woiting the subroutine KPROB. In this subroutine, the 
distance from each atom to the reference (0,0) point in the 
reduced representative impact area was calculated. then all 
atoms at each distance were checked and the atom sputtering 
count was summed. This sum was divided by the total number 


of sputs from all the shots. This quotient was the 
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Conditional probabilty that a given ring at a known distance 
would sputter an atom given the atoms that had already 


Sputtered. 


KPROB performs the above calculations. However, a very 
interesting result is obtained by dividing the conditional 
probabilty obtained above by the number of atoms in a 
Pare@eniart ting. This result was a conditional probability 
that an atom at a given distance from the impact point would 


sputter. 


G. EJECTION ANGLES 


The ejection angle is defined as the angle between the 
atom's velocity vector and the surface normal. A tally was 
kept of each atom's ejection angle, grouped in five dearee 
intervals. This tally was recorded separately for each 
teerope. TO provide a guick analysis of the data, a program 
(airs tO} waS written to print out a histogram of the number 
of sputtered atoms ejected in each interval for each 


isotope. 


Pee TARE ABLE-SHOT SELECTION 


The previous program was only flexible enough to allow 
Perma single=shot run or all 36 shots. As the size of 
the lattice and the ion's energy were increased, the 
computer run time for a set of 36 shots became unmanageable 
ize to hours) . At this point, a variable size shot 
"package'"t was developed. The features included the ability 
to specify the number of shots desired and the starting 


IMpact point of the "package". The final summary for the 


Pu, 





particular package run was printed out. The packages had to 
be combined by hand calculations in order to get the 1e.o)= lll 
information for a particular run. Even so, this procedure 
allowed runs to be made with larger lattices and higher 
energies (i.e. the containment studies) that previously 


would have been prohibitively long. 


meer anDONEY SELECTED IMPACT POINTS 


To obtain more angular dispersion data EOE cae 
low-energy, (100) orientation events, a modification was 
made to the simulation model. Detcoeemnodl ft cation, “Eeplaced 
the systematic method of incrementing the impact point over 
a grid of points in the representative target area with a 
random-shot capability. Again, the IMSL GGUB routine was 
used to generate uniformly distributed, osuedo-random 
numbers. The first random number was labeled the x 
coordinate in the impact area. The second random number was 
labeled the 2z coordinate and tested to determine if it was 
inside the impact area z boundaries determined by the value 
Of the x coordinate. If this test was passed sucessfully, 
these values for the x and 2 coordinates were used. 
However, if the test did not pass (the x and z coordinates 
were outside the impact area), the two original numbers were 
discarded, two more were selected, and the process was 


repeated again. 


DmicommEanadom 2Npact pOint did hot effect the variable - 
shot selaction feature. But, enough random numbers hed to 
be generated to allow for those points that did not fall 


into the impact area. 
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The results of randomly bombarding the impact area will 
be discussed later. This feature was incorporated 59 
response to other investigators suggestion that the grid of 
impact points be compared with a randomly bombarded impact 


ae Cd. 
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IV. RESULTS 


The three lattice orientation models for tne copper 
target were utilized in computer simulation runs in which 
100-eV, 200-eV, 500-eV, and %i1-keV normally incident ions 
bombarded the target. The ion masses varied from 4 anu 
(Het) to over 200 amu (Hg+). The lattice was constructed 
with copper atoms of either the natural mass (63.54 amu) or 
the proper ratio of isotopes (63.93 amu and 64.93 amu). The 
Hattice size varied from a 9x4x9, (100) orientation of 162 


meomses2O a 17K4x17, (100) orientation of 578 atoms. 


Computer runs were made for various permutations of ion 
energy, ion mass, lattic2 orientation, lattice sizes, and 
Betential function. The summary of each run contained the 
humber of atoms sputtered from each impact point, the number 
of times each atom sputtered, a histogram showing the 
angular dispersion of the sputtered atoms for each isotop:, 
the number of sputs for each atom after the folding process, 
the distance from each ring to the impact area reference 
point, the number. of atoms and the sum of tneir sputs in 
SaewmenIng,sand the probability that each ring contained a 


sputtered atom. 


In the sputtering summaries the numbers in the ellipses 
indicate the number of times that particular atom sputtered 
in the simulation. The dark areas, a triangle for the (100) 
anGmwtnli)  OFVentationS and a rectangle for the (110) 
orientation, indicate the location of the target impact 
area. The atom that is struck by the ion is called the 
"DOrimary Knock-on atom" (hereafter PKA). Since different 


seeds were used to generate the random isotope sites, a 
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letter L or H denotes whether the sputtered atom was a 
light or heavy isotope. Wherever possible, the results for 
the Gibson-II and KSE-B potential functions were combined. 
In the figures displaying the sputtering summary and the 
mechanism tracing, some of the lower layers were omitted. 
However, all atoms were included in all calculations. [In 
the 100-eV results all the sputtered atoms are shown. 
Wherever a direction in the lattice is specified, it is 
always with respect to that particular lattice axes. The 


top surface is always an XZ plane. 


pee OX MICROCRYSTALLITE 


Initially, all simulations used a 9x9 lattice Lo 
determine the sputtering ratio as a function of the ion mass 
at 100 eV and 1 keV. The model was modified to include the 
isotope generator, the symmetrical folding of the summary, 
the angular dispersion of sputtered isotopes, and the "ring 


Sputtering" probabilities. 


fe owe Or lentat ion 


The (100) orientation MieCnhoOCcEyVsotalLlate was 
constructed ina lattice consisting of 162 atoms. The 
lattice was 9 atoms in the X and Z directions and 4 atoms 
deep in the ¥Y direction, hence 9x4x9. Figure 3 shows an 
isometric projection of the lattice and its atom numbering 
scheme. Figure 9 shows an expanded view of the grid of 36 


impact points used. 
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a. 100-eV Ions 


Of primary importance was the investigation of 
the angular dispersion of isotopes reported by Wehner [24]. 
His experiments were conducted with 100-eV mercury ions 
bombarding a copper target. These conditions were ideal to 


Simulate with modifications to the existing model. 


The sputtering simulations for a normally 
mieident, 100-EV i10n were run using Het, Net, Art, Cut, 
Mey ket, Aut, and Hgt Tots. Figure 14 shows the 
Seulecerang LEatio aS a function of the ion mass for both 
potential functions. As mentioned in the objectives, a 
decrease in the yield was evident for the heavier ions. The 
KSE-B curve shows a more natural, smooth behavior than that 
Sec oSOn-Lil. 


(1) fot Une Lons. The results when 100-sV 
helium icns impact the target area are shown in figures 15, 
16, and 17. Figure 15 shows the sputtering summaries using 
(a) Gibson-II, resulting in a sputtering ratio of 0.444, and 
Pec on-b, EFEeSUlting in a sputtering ratio of 0.111. The 
predominance of the sputtered atoms in the <100> and <101> 
directions is quite evident. Figure 16 shows he 
Sena'tional probability that a ring of atoms, as a function 
of the ring distance from th2 impact area reference point, 
MenecOontatn dn atom that sputtered for (a) Gibson-II and 
foe KSE-B. The term Loess 0s | sputtering" HEODADLLICY, 
henceforth ring probability, was used as a label to mean the 
above conditional probability. The relative amplitudes of 
the ring probabilties could be misleading becaus¢ the 
calulation does not take into effect the number of atoms in 
each ring. Again, the predominance of rings that contain 


meetseinme the <100> and <101> directions is evident. The 
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peak of the KSE-3B curve was closer to the center than the 
Gibson-II peaks, which were spread out all the way to the 
lattice edges. Figure 17 shows the angular dispersion 
results for (a) Gibson-II and (b) <xKS&-8. The particular 
seed used in the psuedo-random number generator to locate 
the isotope random sites was such that only light isotopes 


sputtered in this case. 


(2) Neon Ions. The results when 100-e€V neon 
1ons impact the target area are shown in figures 18, 19, and 
20. Figure 19 shows that for the light mass of the neon 
lon, atoms on the lattice boundary are sputtered using (a) 
Gibson-II and (b) KSE-B. Gibson-II resulted in a sputtering 
meelOoewroQrt 0.8989, while the sputtering ratio for KSE-B was 
0.806. Figure 19 shows that the Net sputtered atoms were 
mmitnmer aWdy from the lattice center than the Het, for both 
feyeectpbson-Ii and (b) KSE-B. This was confirmed by the 
increased number of next-nearest neighbor atoms (in the 
<101> direction) sputtering. Figure 20 indicates that some 
of the heavier isotopes sputtered for (a) Gibson-II and (b) 
KSE-B. There were not enough sputtered atoms to notice any 


trend in the isotore angular dispersion. 


(3) Argon Ions. The argon Simulations were 
Heed tO ianvestigate the "reflective collision" theory. 
Argon ions of 100-eV energy bombarded the target area using 
both of the potential functions. In the random isotope site 
locations, the PKA (atom 22) was a heavy isotops, while 
feee (52, 58, and 67) of the four atoms directly underneath 
were light isotopes. Four simulations were run using ¢ach 
Soe tn] §potential functions; (1) the random isotope sites, 
(2) forcing the PKA to a light isotope, (3) forcing the PKA 
*<O a heavy isotope and all four atoms underneath to a light 
moomos=, ana (4) forcing the PKA to a light isotope and all 
four atoms underneath to a heavy isotope. The remaindar of 


the lattice isotopes were unchanged. 
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Figure 21 shows the sputtering summari2s 
for all four combinations using (a) Gibson-II and (b) KSE-B. 
In all four cases, identical atoms sputtered independent of 
isotope locations. Gibson-II resulted in sputtering ratios 
See ts0>0, While KSE-B resulted in sputtering ratios of 
Meee. Figure 22 illustrates the ring probabilty in all 
four cases for (a) Gibson-II and (b) KSE-B. Plame 22 
lllustrates the angular dispersion of ejected atoms, using 
Seesen-i1f, for (a) case (1), (b) case (2), (c)”case (3), and 
(d) case (4). No light isotopes were sputtered. Pigure 24u 


illustrates the same results for KSE-B. 


Since at least four atoms were ejected into 
a different five degree window for Gibson-II, the sputtering 
mechanisms were checked. The most active event occurred 
when the impact area was struck by the ion at impact point 
Pere i imeact displaced 0.40 LU in the X direction and 0.60 
fueeeen ene 27 direction from the center of the PKA, hence 
4060). Four atoms wer? sputtered at this POUnt A 
Ssingle=impact program was run at this point with the output 
[Mietynig tne positions, velocities, kinetic energy, and 


potential energy of each atom with a total energy greater 


man 0.1 eV. The listing was printed once ever five 
timesteps. A timestep is the time interval between dynamic 
calculations. It is variable and determined as the time 


required for the fastest moving atom to move 0.1 LU from its 
previous position. In real time the timestep is on the 
order Or 10-15 seconds. Pigure 25 illustrates aton 
Meajececories in the sputtering of atoms 17, 32, 57, and 7/7. 
ipGases §(3) and (4), the sputtering mechanisms were the 
same displaying only a slight difference in the ejection 


angles and the velocities of the sputtered atoms. 
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The argon ion collides with atom 22 # £(PKA) 
and atom 27 by timestep 5, pushing them down and under atoms 
17 and 32 respectively. The ion continued downward and 
collided with atom 67 by timestep 15. The ion reflected 
back up through the hole vacated by atoms 22 and 27. By 
timestep 15, atoms 17 and 32 were moving upward. Atom 67 
moved down colliding with atoms 107 and 112 by timestep 29. 
Atoms 107 and 112 squeezed [21] under atoms 102 and 117 
respectively. Atoms 102 and 117 were moving up by timestep 
30, and struck atoms 57 and 77 from below by timestep 40. 
By this time, atoms 17 and 32 had moved far enough above the 
lattice with sufficient upward velocity to be classified as 
eapucs'. Atom 17 was ejected at an angle of 33.56° to the 
surface normal with an upward velocity of 2924 m/sec for 
ersen(s), and at an angle of 33.369 with a velocity of 2974 
m/sec for case (4). Atom 32 was ejected at an angle of 
Eemono° With an upward velocity of 2926 m/sec for case (3), 
monde at an angle of 33.569 with a velocity of 2918 m/sec for 
case (4). Atoms 57 and 77 continued upward through the 
holes left by the sputtered atoms 17 and 32. As atom 57 
passed tnrough this hole its trajectory was bent closer to 
the surface normal by atoms 12, 16, and 21. Likewise, atom 
tewas crefilected upward by atoms 28, 33, and 3/7. BY 
timestep 75, atoms 57 and 77 were far enough above the 
lattice surface with sufficient upward velocity to ope 
Classified as "sputs". Atom 57 was ejected at an angle of 
18.129 with an upward velocity of 2967 m/sec for case (3), 
and at an angle of 17.75° with a velocity of 2882 m/sec [or 
case (4). Atom 77 sputtered at an angle of 25.049 with an 
upward velocity of 2799 m/sec for case (3), and at an angle 
of 24.759 with a velocity of 2859 m/sec for case (4). 
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Tt is noted that when the PKA was a light 
isotope, all of the sputtered atoms were ejected closer to 
the surface normal. Also noted is the remarkable symmetry 
of ejection angles and velocities of the sputtered atoms. 


(4) Copper [Ions The results when 100-eV 


copper ions impact the target area are shown in figures 26, 
a, and = 28. Figure 26 shows the sputtering summaries for 
(a) Gibson-II, resulting in a sputtering ratio 0.944, and 
foeeeaSh-B, resulting in a sputtering ratio of 1.111. 
Gibson-II Simulations with Art and Cut were the only runs 
that had atoms sputter from the second layer. Also, these 
two ions resulted in the highest sputtering ratios for both 
potential functions. These results further indicate a 
Maximum sputtering ratio occurs near an ion/atom masS fratio 
Of unity. Figure 27 illustrates a continued spread outward 
See tene Ling probability for (a) Gipson ian cus D)peeno 5-3. 
Figure 28 illustrates the angular dispersion results for (a) 
G@bson-TI and (b) KSE-B. 


(5) ihe ypcon LOons. The results when 100-eV 
krypton 10ns impact the target area are shown in figures 29, 
Oo, and 31 Figure 29 shows the sputtering summaries for 
fey Grbpson-fi, resulting in a sputtering ratio of 0.694, and 
[jeeetoe-B, resulting in a sputtering ratio of 1.139. dere, 
a Maximum in the sputtering ratio was reached for XSE-B, 
while the sputtering Pa to Ons Gibson-II dropped 
@mastiecaltiy. Figure 30 illustrates a continued outward 
Spread of the ring probabilities for (a) Gibson-II and (b) 
KSB=-B. ft is noticed that the outward spread still occured, 
even though the Gibson-II yield decreased. Figure 31 
illustrates the angular dispersion results for (a) Gibson-II 
and (b) KSE-B. 
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(6) Xenon Tons. The results when 100-9V xenon 
ions impact the target area are shown in figures 32, 33, and 
34. Figure 32 shows the sputtering summaries for (a) 
eeeecon-ii, resulting in 2 sputtering ratio of 0.556, and (b) 
KSE-B, resulting in a sputtering ratio of 1.028. Figure 33 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. The outward spread of the ring probability peaks is 
still seen for KSE-B, while Gibson-II remained constant. 
Figure 34 illustrates the angular dispersion results for (a) 


SPoson-It and (b) KSE-B. 


(7) Gold 


ions impact the target area are shown in figures 35, 36, and 


ons. The results when 100-eV gold 


37. Figure 35 shows the sputtering summaries for (a) 
Seeoscon-li, resulting in a sputtering ratio of 0.528, and (b) 
Rop-b, resulting in a sputtering ratio of 0.778. Figure 36 
illustrates the ring probabilities for (a) Gibson-II and (b) 
KSE-B. Gibson-II has settled out with approximately equal 
probability of sputtering from atoms in the first two rings 
fmieche <I01> direction. In the case of KSE-B, a shift fron 
the first and second rings in the <100> direction toward the 
secend ring in the <101> direction is noticed. Figure 37 
illustrates the angular dispersion results for (a) Gibson-II 
and (b) KSE-B. 


(8) Mercury Ions. These Hgt bombardment 
results are compiled from the sets of Simulation runs’ used 
to investigate the "reflective collision" theory. Figure 38 
shows the sputtering summaries for all four cases, in which 
the masses of the PKA and the four atoms underneath are 
Varied. The results for (a) Gibson-II were identical in all 
meueeecases, vcesulting in a sputtering ratio of 0.528. The 
Sottreting EFcatio was 0.743 using KSE-B in a Lattnce 
Genmscleting of atoms with only the natural mass. fhe results 


for (b) KSE-8 depended on the order of the heavy over light 
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interaction. When the PKA was a heavy isotope, shown in 
figure 38(b), the sputtering ratio was 0.7229 However, when 
the PKA was a light isotope, atom 30 sputtered an additional 
time (denoted by A) and the sputtering ratio became 0.750. 
Figure 39 oiusera cess theesring probabilties for -~ (a) 
Gibson-II, which remained relatively unchanged, and (b) 


foe, which still displayed a shift to the second ring in 


the 101 direction. Figure 40 illustrates the angular 
dispersion results for Gibson-II. Figure 41 shows the 
aiauear dispersion results for KSE-B. Again, atom 30 


Sputtered an additional time. The additionai ejection angle 
is denoted by A. 


om 1-keV Results 


The 1i-keV events were obviously not contained in 
the 9x4x9 microcrystallite. However, the simulations were 
run to compare the sputtering ratios with exverimental 
results. Also, the dependence of the sputtering on ion mass 
Was investigated. Useful information can be obtained about 
the ring probabilties and isotope angular dipersion at this 
higher energy. The 162 atom, 9x4x9 microcrystallite, with 
the correct ratio of copper isotopes at random sites, was 
used. Even though all atoms were used in the calculations 
and the printout provided information on all atoms, only the 
top layer results are shown in the figures. pve (eee te 
worse case (Ciy7 AG GyeGibSon—-iLl), more than 86 per cent of 


the sputtered atoms were from the top layer. 


The Seer 21g Simulations Bone normally 
mietmdent, Im=keY 10nNS were run using Het, Net, Art, Cut, 
Kr+, Xet, Aut, and Hg+ ions. Figure 42 shows the sputtering 
Boeowasea L£UNCLIOnN Of the ion mass for both potential 


Bunce LOnS. As expected, the sputtering ratios have 
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increased. However, they are no longer smooth and 
continuous, but exhibit a randomness, as if each was 


independent of the others. 


All the 1-keV simulations used the same random 
isotope sites. The letters L and H, denoting light or heavy 
isotopes, appear only in the Cu/Net, KSE-B summary (more 
different atoms sputtered - but not the highest sputtering 


ma LO) . 


(1) Helium Ions. The results when 1-keV heliun 
ions impact the target area are Shown in figures 43, 44, and 
45, Figure 43 shows the sputtering summaries for (a) 
Seeson-ii, resulting in a sputtering ratio of 1.914, and (b) 
fee resulting in a sputtering ratio of 2.833. ite oot h 
cases, boundary atoms sputtered. Figure 44 illustrates the 
memgeOrODabilities £or (a) Gibson-IIT and (b) KSE-B. 8Soth 
display a large action-reaction effect, as if a bomb were 
dropped at the center and everything out from the center was 
blasted upward. Figure 45 lllustrates the angular 
dispersion results for (a) Gibson-II and (b) KSE-B. 

(2) Neon Ions. The results when 1-keV neon ions 
impact the target area are shown in figures 46, 47, and 48. 
Figure 46 shows the sputtering summaries for (a) Gibson-II, 
me=sule ing In a Sputtering ratio of 3.686, and  (b) KSE-B, 
Zometing 2n a Sputtering ratio of 4.500. Host of the 
Surface atoms were sputtered, indicating that the energy /( 


Or momentum ) spread out over the entire lattice surface. 


Figure 47 illustrates the ring propabilities for (A) 
SE soson-r£t and {(b) KSE-B. The KSE-B simulation indicate2s a 
Spreading outward of the sputtered atoms. The Gibson-II 


Simulation not only exhibits increased sputtering of atoms 
further from the center, but it also clearly shows a snift 
in the number of sputtered atoms from those at 2.83 LU in 
mmem—eror> direction to those at 3.16 LU in the <103> 
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ieee eectlon. AS mentioned earlier, this could be artificial 
because the ring at 3.16 LU is the first ring to contain 
eight atoms. Figure 48 illustrates the angular dispersion 
results for (a) Gibson-II and (b) KSE-B. 


(3) Argon Ions. The results when 1-keV argon 
ions impact the target area are shown in figures 49, 50, and 
Bik Figure 49 shows the sputtering summaries for (a) 
measson- hh resulting in a Sputtering ratio of 3.806, and ({(b) 
fees, cesulting in a sputtering ratio of 4.611. Pigure 50 
Piewwstrates the ring probabilities for (a) Gibson-II and (%) 
PoE-B. The Gibson-II simulation continued to increase in 
the outward direction at the expense of the atoms at 3.16 LU 
fae the <103> direction. However, the KSE~-3B simulation 
Meversed 1tS trend and the ring probability for the close-in 
rings (1.41 LU and 2.00 LU) increased at the expense of the 
mengeat 2.63 LU andsthe rings further out from ‘the center. 
Figure 51 illustrates the angular dispersion results for {a) 
PaosOn-Lf and KSE-3. 


(4) Copper Ions. The results when 1-keV copper 
10nNS impact the target area are shown in figures 52, 53, and 
54. Figure 52 shows the sputtering summaries for (a) 
Peoeon-tt, resulting in a Sputtering ratio of 3.028, and (b) 
moe-oc, cCesulting in a sputtering ratio of 3.972. These low 
Sputtering Catios not only disagree with a monatonic 
increase in the yield with mass, but are drastically low for 
the maximum to occur at a ion/fatom mass ratio of unity. 
Figure 53 illustrates the ring probabilities for (a) 


Gibson-II and (b) KSE-B. The Gibson-II simulation showed an 


ubexpected change. The ring probability increased for the 
mergceacme.O0 LU (<i00> direction) and 3.16 LU (<103> 
g@irect ion) . This trend was also displayed in the KSE-B 
Samnelacton, Where it starts to occur in the OS ees 
Simulation. Figure 54 illustrates the angular dispersion 


Beovetse=or (a) Gibson-IE and (b) KSE-B. 
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(5) Krypton Ions. The results when 1I-keV 
Krypton ions impact the target area are shown in figures 55, 
Beye and 5/7. Figure 55 shows the sputtering summaries for 
foyeecabson-ii, resulting in a sputtering ratio of 2.833, and 
fee sSe-B, resulting in a sputtering ratio of 3.917. Note 
maae the Sputtering ratio, for both potential functions, 
continues to decrease with increased ion mass. Figure 56 
illustrates the ring probabilities for (a) Gibson-II and (b) 
hRoa-B. Both potential function Simulations continue to 
exhibit the trend to establish a peak at 2.00 LU in the 
MimOeeratrection and at 3.16 LU in the <103> direction. 
However, the KSE-B simulation shows significantly larger 
meng probability at 3.16 UJ. Figure 57 illustrates the 


Bmoulbat dispersion results for (a) Gibson-II and (b) KSE-B. 


(6) Xenon Ions. The results when 1-keV xenon 
mons impact the target area are shown in figures 58, 59, and 
0. Figure 58 shows the sputtering summaries for (2a) 
[eescOn-il, Eesulting in a sputtering ratio of 2.333, and (b) 
meno, resulting in a sputtering ratio of 3.472. Both 
potentiai functions continued to show a decrease in the 
Sputtering ratio with an increase in the mass of the ion. 
Figure 59 illustrates the ring probabilties fOr (a) 
Seoson—-iE and {b) KSE-B. All the previous trends, except an 
increased ring probability at 3.16 LU were shattered. 
Figure 60 illustrates the angular dispersion results for (a) 


Gibson-ITI and (b) KSE-B. 


‘7p Gold ions 


1ons impact the target area are shown in figures 61, 62, and 


The results when 1-keV gold 


Os. Figure 61 shows the sputtering summaries for (4) 
Gibson-II, resulting ina Sputtering ratio of 2.778, and (b) 
KSE-B, resulting in a sputtering ratio of 4.167. A definite 
increase in the sputtering ratio occured. Figur: 62 


illustrates the ring probabilities for (a) Gibson-II and (b) 
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KSE-B. No new trends were noted, except that the fring at 
B10 LU continues to dominate and the other fring 
probabilities fluctuate. Figure 63 illustrates the angular 


dispersion results for (a) Gibson-II and (b) KSE-8. 


P (8) NemeMEy fons. The results when 1-kev 
mercury ions impact the target area are shown in figures 64, 
wopand 66. Figure 64 shows the sputtering summaries for 
fpeeGibson-ii, resulting in a sputtering ratio of 2.778, and 
asyekoo-B, resulting in a sputtering ratio of 4.278. These 
Simulations resulted in no change in the sputtering ratio 
PieGriOseon-if, and only a slight increase for «SE-B. Figure 
65 illustrates the ring probabilities for (a) Gibson-II and 
(b) KSE-B. Only insignificant changes occured compared to 
meme gold icn results. Figure 66 illustrates the angular 


dispersion results for (a) Gibson-II and (b) KSE-B. 


on 2-keV Results 


The 2-keV ion simulations were the earliest runs 
made. No isotope information waS obdtained. AS a 
consequence, the only results are the sputtering ratio 
versus ion mass curve, figure 67, for both potential 
functions. They exhibit an even more erratic behavior than 
did the 1-keV simulations. The sputtering summaries 
frequently show atoms sputtered from the third layer, and 
rarely an atom sputtered from the fourth layer. At certain 
impact points aS many as 25 atoms sputtered. The small, 
9x4x9 microcrystallite was effectively plasted apart by the 


2-keV ions. 
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Zee 0) Orientation 100 eV 


tLReweonroru cevoted "te. (110) orientation simulation 
runs was very limited due to the higher priority assigned to 
feme (100) Orientation runs. As a result, only the argon and 
mercury ion Simulations at 100 eV were run. Still needing 
investigation are the sputtering ratio as a function of ion 


mass and sputtering mechanism traces. 


The (110) orientation NMECFOGE YS tall ite was 
constructed ina lattice consisting of 123 copper isotopes. 
In the random isotope site locations, the PKA (14) was a 
imght  2Setope, as were three (33, 36, and 37) of the four 
atoms underneath the PKA. Figure 4 shows an isometric 
projection of the lattice and its numbering scheme. Figure 
10 shows an expanded view of the grid of 36 ‘impact points 
used. A 9x6x9 lattice was used with no increase in computer 


time required to complete a run. 


Two sets of simulations were run to investigate the 


"reflective collision" theory. Four runs were mad@ in each 
set of Simulations for each potential function. The impact 
area was bombarded by 100-eV, argon and mercury ions. five 


four runs consisted of; (1) the random isotope sites (PKA 
was a light isotope), (2) forcing the PKA to ae heavy 
Poot ope, (5) f£oOrcing the PKA to a light isotope and all four 
atoms underneath to a heavy isotope, and (4) forcing the PKA 
to a heavy isotope and all four atoms underneath to a light 
isotope. The remainder of the lattice isotopes were 


unchanged. 
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an AEGON Lons 


The results when 100-eV argon ions impact the 
target area are Shown in figures 68-75. Figure 68 shows the 
Sputtering summaries using Gibson-II for all four cases. 
Case (1) and case (3), shown in figure 68(a), resulted in a 
PeeiarerLing ratio of 0.38389. Case (2) and case (4), shown in 
meouee 66 (b), resulted in a sputtering ratio of 0.861. A 
predominance of sputtered atoms in the <100> and <101> 
Mecections is evident. PugGise moo metiliustmates, “the | ring 
prebabilities for (a) case (1) and case (3), and (b) case 
(2) and case (4) uSing Gibson-II. Notice that an atom at 
the boundary (in the <101> direction ) svuttered. Figure 
70 illustrates the angular dispersion results for (a) case 
fie and (b) case (3) uSing Gibson-II. Figure 71 shows the 


Same information for (a) case (2), and (b) case (4). 


Figure 72 shows the sputtering summaries using 
KSE-B for all four cases. Case (3) resulted in a sputtering 
ratio of 0.750, while case (1) sputtered atom 26 one 
additional time (denoted by A) and had a sputtering ratio 
Of 0.778. Case (1) and case (3) are shown in figure 72 (a). 
Case (2) resulted in a sputtering ratio of 0.778, while case 
(4) Sputtered atom 8 one additional time and had a 
Sputtering ratio of 0.805. Case (2) and case (4) are shown 
moe f£igure 72(b). Again, it is evident that only atoms in 
the <100> and <110> directions sputtered. Figure us 
MmmeistGates the ring probabilities, using KSE-B, for (a) 
case (1) and case (3), and (b) case (2) and case (4). The 
close-in trend, Similar to-that of Gibson-II, is evident, as 
are the sputtered atoms at the boundary in the <i101> 


dEectlon. Figure 74 illustrates the angular dispersion 
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meotlbes, USINg KSE-B, for (a) case {1) , and (b) case (3). 
Figure 75 shows the same information for (a) case (2), and 
(b) case (4). 


b. Mercury Ions 


The results when 100-eV mercury ions impact the 
target area are shown in figure 76, 77, and 78. Figure 
76(a) shows the sputtering summaries for all four cases, 
meng  GibSoOn-ii. Case (1), case (2), and case (4) resulted 
Mims puttering ratio of 0.278. Case (3) sputtered atom 13 
one additional time (denoted by A) and had a sputtering 
maeeo Of 3.05. Figure 76(b) shows the sputtering summaries 
for all fOwrEuN cases, Using  KSE-B, which resulted in 
PMiecering ratios of 0.222. Figure 77 illustrates the fring 
BEODaDilitcies for (a) Gibson-II and (b) KSE-B. Figure 78 
illustrates the angular dispersion results for (a) Gibson-II 
ea (Db) KSE-B. 


No significant trends were noted in the (110) 
orientation. However, the ease with which atoms sputtered 


mim the <10l> direction is evident. 


Be (111) Orientation 100eV 


Like the (110) orientation, thea (111) orientation 
did not receive extensive attention. Both the sputtering 
ratio as a function of the ion mass, and sputtering 


mechanism traces need further investigation. 


The (ol tal) orientation microcrystallite is 
memosenueted in a 9x4x9 lattice consisting of 163 atoms. [In 
the random isotope site locations, the PKA (22) was a light 


Isotope, as were three (58, 63, and 67) or the rour atoms 
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underneath the PKA. Figure 5 shows an isometric projection 
of the lattice and its numbering scheme. Figure 11 shows an 


expanded view of the grid of 21 impact points used. 


Two sets of Simulations were run to investigate the 
Meeftlective collision" theory. Four runs were mad2 in each 
eee Of Simulations for each potential function. The target 
was bombarded by 100-eV argon and mercury ions. Tie Otc 
runs consisted of; (1) the random isotope sites (PKA was a 
light isotope), (2) forcing the PKA to a heavy isotope, (3) 
forcing the PKA to a light isotope and all four atoms 
underneath to a heavy isotope, and (4) forcing the PKA to a 
Meavyy isotope and all four atoms underneath to a light 
isotope. The remainder of the lattice isotopes qere 


unchanged. 


am PArgoOn Lons 


The results when 100-eV argon ions impact the 
mifeget area are Shown in figures 79-83. Figure 79 shows the 
sputtering Summaries HO ta GlUsOn- Phe resulting in 
Perec ering Eatios of 0.857, and (b) KSE-B, resulting in 
Semccer ing ratios of Q0./62. For KSE-B, case (2), atom 16 
sputtered one additional time (denoted by A), resulting in 
MEemecering ratio Of 0.3810. Figure 80 illustrates the ring 
Peepabilities £0r (a) Gibson-II and (b) KSE-B. It was 
unexpected that, for KSE-B, all of the sputtered atoms are 
iameme < l01> direction in the ring at 2.83 LU. Figure 81 
illustrates the Gibson-II angular dispersion results for (a) 
case (1) and case (3), and (b) case (2) and case (4). 
Figure 82 illustrates the KSE-B angular dispersion results 
for (a) case (1) and (b) case (3). Figure 83 shows the same 


results for (a) case (2), and (b) case (4). 
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Db. ) Hencury [ons 


The results when 100-eV mercury ions impact the 
target area are shown in figures 84-89. Figure 84(a) shows 
the sputtering summaries, using Gibson-II, for case (1), 
case (2), and case (4), which resulted in sputtering rcatios 
Omee0.476. Case (3) had a sputtering ratio of 0.381 and iis 
shown in figure 84 (b). Figure 85 shows the sputtering 
summaries using Pniewe ho —Ei. Case (3) resulted in a 
Sputtering ratio of 0.762, while case (1) sputtered two 
additional atoms (12 and 14 - denoted by A) and has a 
BmECeCEINng Fatio of 0.857. Case (1) and case (3) are shovwn 
Mmieerigure 85(a). Case (2) and case (4), shown in figure 
eeepoy, Nave Sputtering ratios of Cmts 2. Figure 86 
mis trates the ring probabilities for (a) Gibson-II and (5) 
KSE-B. Again, it is evident that the only sputtered atoms 
mae th the ring at 2.83 LU, in the <101> direction. Figure 
87 illustrates the Gibson-II angular dispersion results for 
(a) case (1), case (2), and case (4), and (b) case (3). 
Figure 88 illustrates the KSE-B angular dispersion results 
for (a) case (1) and (pb) Snes (3). Figure 89 shows the same 


information for (a) case (2), and (b) case (4). 


fee GARGER LATTICES 


The erratic behavior of the sputtering ratio asa 
mine lion Or the On mass, indicates that the 
microcrystallite is not containing the event. Originally, 
tcontainment" was defined as "no boundary atoms sputtering." 
Based on this definition, none of the 100-eV events vere 
@eomrained in the (100) orientation or (110) orientation 9x9 
Meeroecrystallites. However, in the (111) orientation, 


100-eV events were contained. The sputtering summaries for 
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mercury tons pombarding a  t3X6x13, CUTOy=— Deaen tation 
MmerOctyStallite containing 255 atoms are shown in figure 
90. Figure 90(a) is for Gibson-II, while figure 90(b) is 
mom KSE=B. Comparison of figure 90 with figure 1% 
Seontirmed that the 9x5x9 microcrystallite did not contain 
the 100-eV, mercury event. Based on the original definition 
@emcontainment, the 13x6X13 microcrystallite did not contain 


the Gibson-II Simulation run. 


In the (111) orientation, the sputtering summari2s for 
mercury 1ons bombarding a 13x4x13 microcrystallite 
Someaining 339 atoms are shown in figure 91. FPigure 91 (a) 
mom fOr Gibson-li, while figure 91{b) is ZO Koad 
Momparison of figure 91 with figures 84 and 85, indicate 
that even though the 9x4x9 lattice had no boundary atoms 
Sputtered, the sputtering ratio increases when a 13x4x13 


lattice is used. 


1. 100-eV¥ Containment 


Based on the above information, the definition for 
"containment" was changed to: "an event is contained when an 
increase in the mMicrocrystallite size does not result in any 
further increase i sene — yterd. Barlier simulations 
Pweeoea an darger yields for KSE-B. For this reason, KXS5E-B 
was used to investigate the size microcrystallite needed to 


contain a 100-eV event. 


The sputtering ratio as a function of the ion mass 
Was examined for the (100) orientation. Figure 92 shows the 
meats for Ix4x9, tix&xil, and 11x6x11 microcrystallites. 
Additionally, simulations were run for copper (largest 
Se@ereLing EFatio) and metccury ions for 13x4x13, 13x6x13, and 
17x4x17 microcrystallites. Based on the new definition of 


containment, 100-eV¥V helium and xenon events are contained in 
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aeeoxk4 xo Microcrystallite because an 11x4x11 microcrystallite 
produces the same yield. However, the argon, copper, and 
krypton events were not contained. Further increase of the 
microcrystallite size demonstrated that the 100-eY, copper 
event is contained by a 13x4x13 microcrystallite. Figure 92 
also illustrates that the maximum sputtering ratio cccurs 
when the lon/fatom mass ratio is approximately unity. Figure 
fee Shows the same information in a plot of sputtering ratio 


versus lattice size for all eight ions. 


Zee ss X17 Mrcrocrystallits 


Asuming that all 100-eV events are contained by a 
17x4x17 microcrystallite, the effect of increasing the ion 
energy was investigated. The reason for this investigation 
was to determine how energetic an event could be contained 
by a 17x4x17 microcrystallite. And to compare higher energy 
events With experimental results. The reason for not 
increasing the lattice size without limit was obviously to 
keep the computer run time within reason. Again, only KSE-B 


with mercury and copper ions were investigated. 


a. Copper Ions 


The results when 100-eV copper ions impact the 
target area are Shown in figures 94, 95, and 96. Figure 94 
Shows the sputtering summary with a sputtering ratio of 
ieeoo9 . Containment for the 100-eV copper event had already 
been established. Figure 2S) illustrates the "atom 
sputtering PRepaa lity”, henceforth called the 
Seuedanittaty™. The "probability" is the "ring sputtering 
probability" divided by the number of atoms in each ring. 
The dimensions for "probability" is sputs per atom divided 


by the total sputs at a fixed distance from the impact area 
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reference point. Figure 95 indicates an enhancement of 
sputtered atoms in the <100> and <101> directions. Figure 96 


illustrates the angular dispersion results. 


Results for the 200-eV copper ion simulations 
are shown in figures 97, [6 and - 99. Price 97 jee ne 
sputtering summary with a sputtering ratio of 2.8389, shows 
that boundary atoms are sputtered, and that the event 
@efinitely is not contained. Figure 98 illustrates the 
Mprebability"™ results and indicates that an atom farther out 
had increased chances of sputtering. Figure 99 illustrates 
the angular dispersion results. A bias toward IE ale alae 


isotopes closer to the normal is evident. 


The results (ONG the 500-eV copper 1 
Simulations, with a sputtering ratio of 4.028, are shown in 
mrqures 100, 101, and 102. Again, boundary atoms are 
sputtered and a slight shift outward of the "probability" is 
noticed. Figure 102 illustrates the angular dispersion 
results which showed that more light isotopes are ejected 


more nearly normal to the lattice surface. 


the 1-keV, copper ion results with a sputtering 
Bacio Of 7.083, are shown in figures 103, 104, and 4105. 
Figure 103, the sputtering summary, shows a very large 
number of the surface atoms sputtered. Boundary atoms in 
every direction are sputtered. This increased spréading 
Outward was confirmed in figure 104, which shows the 
"probabilities" relatively level. Puce | Wo tl iustraces 


the angular dispersion results. 
Gavaousiy, the 200-ev, S00-eV, and I-keVY cooper 


events are not contained, even 138 a 1x7 


Meccocrystallite. 
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b. Mercury Ions 


The results EOT the TOU =Sy SS nereury, -1on 
Simulations are shown in figures 106, 107, and 108. Again, 
the 11x4x11 lattice has already been determined to contain 
the 100-eV, mercury event. Pigure 106 shows the sputtering 
summary with a sputtering ratio of 0.833. Ea Gu bom 
miaustrates the "probability", and like the copper 
Simulation, was predominant in the <101> direction. Figure 


108 illustrates the angular dispersion results. 


The 200-eV mercury Lela results, with a 
meer ering ratio of 1.722, are shown in figures 109, 110, 
eed tt 1. Figure 109, the sputtering summary, shows an 
outward spreading of the sputtered atoms. However, no 
boundary atoms were sputtered and the 200-eV, mercury event 
is tentatively classified as "contained". Pagiure 110 
Semitems the Outward shift in "probability". Figure 11% 
illustrates the angular dispersion results, which indicate a 
greater number of heavy isotopes sputtering than the isotope 


Mati1O Would justify. 


The 500-eV mercury Lem EeSstilts jan tend 


Per ering Fatio Of 3.667, are shown in figures 112, 113, 


and 114. Figure 112, the sputtering summary, indicates that 
boundary atoms are sputtered. PFigure 113 illustrates the 
continued outward spreading of sputtered atoms. Figure 114 


illustrates the angular dispersion results, and the number 


of heavy to light isotopes closer to the correct ratio. 


The 1-keV mercury ion results, with a sputtering 
Berommote 6.139, are shown in figures 115, 115, and 117. 
Figure 115 shows that boundary atoms are sputtered. figure 
iImoeshows the "probability" still spreading outward. Figure 


117 illustrates the angular dispersion results. 
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MicwerpEObability™ ey tigures for both the copper 
and mercury ions are remarkably similar and display the same 


trends. 


c. Sputtering Mechanism Traces 


The containment studies indicated that, as the 
lattice size was increased, the sputtering ratio increased. 
Investigation into the similarities of the 17x4x17 and 9x4&x3S 
microcrystallites, indicated that atom 108 Cl xGcaiy ) 
Sputtered twice in 36 events. This corresponds to atom 40 


mien] 9s4x9 microcrystallite, which did not sputter in 36 


events. Sputtering mechanism trace runs were performed to 
determine the sputtering mechanism involved in the 
Sputtering of atom 108. Atom 108 sputtered when ions 


mmipacted at points 50 ( impact point displaced 9.50 LU in 
fie 2 direction only) and 60 ( impact point displaced 0.6 LU 


maetne 2 direction only). 


Figure 118 illustrates the atom trajectories in 
me] Sputcering of atom 108 from impact point 50 ( the 
mechanism is the same for impact point 60). The copper ion 
Struck atoms 74 (PKA) and 91 forcing them down by timeste 
15. AS atom 91 was moving downward, atoms 99, 100, and 108 
prevented atom 91 flying off in the positive Z direction, 
and channeled it downward. By timestep 25, atom 91 had 
collided with atom 224, reversing its direction. At 
timestep 30, atom 91 was forcing atom 108 upward from 
underneath. Atom 91 transfered almost all of its energy 
(momentum) to atom 108 and had less than 0.1 eV by timeste 
45. As atom 108 was going upward, it also had a velocity in 
the positive Z direction of greater than 6000 m/sec. Atoms 
116, 117, and 125 bent atom 108 upward quite radically. The 
Memeecity inh the Z direction had been reduced to less than 


700 m/sec by timestep 50. By timest2p 60, atom 108 was far 
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enough above the lattice surface, with sufficient velocity 
to be classified as a sput. Atom 108 was ejected at an 
emote of 19.239 to the surface normal with an upward 


velocity of 3206 m/sec. 


Quite evident is the role atoms 116, 117, and 
125 have in channeling atom 108 upward. In the 9xticg9 
lattice, atom 4O ( a boundary atom) corresponds to atom 108. 
Mma nothing equivalent exists for atoms 116, 117, and 125. 
Investigation was then made to determine whether an atom 
corresponding to atom 103 could be sputtered in the 9x4x9 
lattice by using the same impact point on a different PKA. 
[eetnts Simulation, atom 31, corresponding to atom 108 in 
the larger lattice, was expected to sputter. Figure 119 
illustrates the atom trajectories when the copper ion struck 
the new PKA (atom 13) at impact point 50. Atom 31 was 
ejected at an angle of 20.179 with an upward velocity of 
3127 m/sec. The mechanism is the same, and the importance of 


boundary atoms is very evident. 


Additional runs were made on the 9x4x9 lattice 
Bom barding atom 13 at impact points 50 and 060. At both 
points, comparison is made between runs witn (1) only atoms 
of the natural mass, , and natural mass, and (2) PORG ins 
atom 22 (91) to a heavy isotope and atom 67 (224) to a light 
isotope. In all four cases, the sputtering mechanism 1s tne 
same. The differences in the ejection angles and upward 
Mebocities, for each impact point, are less than 0.4° and 75 


m/sec respectively. 


3. Jon/Atom Mass 


Sigmund {28] proposed a theory that the maxinun 
Sputtering vield increases as a function of energy. 


However, results of this investigation strongly indicate a 
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dependence on the ratio of the ion mass to the target atom 
mass. A rearrangement of Sigmund's theoretical curves 
(figure 6b, ref. 28) to show the sputtering ratio as a 
ire tion Of the ion/atom mass ratio results in figure 120. 
The curves shown by Sigmund are predicted for 209-eV, 
600-eV, 1-keV, and 5-keV ions bombarding an aluminum target. 
Figure 120 illustrates that Sigmund's theory would predict 
mime Sputtering ratio to increase as a function of the 
ion/fatom mass ratio. Figure 121 shows experimental results 
obtained by Wehner and is presented by Sigmund ( figure 12, 
mi, 26 ) to support his theory. This figure has. been 
replotted to show the sputtering ratio as a function of the 
Hon/atom mass ratio. The results were obtained  d5y 
bombarding a copper target with 100-eV, 200-eV, and 600-eV 
ions. Experimental results plotted in this manner show a 
decrease in the sputtering ratio with an ion/atom mass ratio 
greater than unity, in disagreement with the theory proposed 
by Sigmund. Results obtained in the present investiga*ion 
completely agree with the shape of experimental data shown 


mi c£igur= 121. 


See kANDOM IMPACT POINTS 


In response to other investigators! suggestion that the 
target impact area be bombarded randomly rather than the 
PyeccMatic grid of impact points, the computer program was 
Slightly modified. Since additional data were desired for 
mercury events, Hgt was used as the incident ion. The 
KSE-B potential function was used. #arlier results reveaied 
that the 100-eV, mercury ion event waS contained in an 


11x4x11 microcrystallite. 





Figure 122 illustrates the (100) orientation impact area 
that was struck at the 144 randomly selacted impact points. 
A few gaps are evident, but the randomness iS very apparent. 
Pegures 123-128 show a comparison of the results from a 
randomly tombarded impact area and the results from the grid 
of impact points. Since the number of events were not the 
Same, the sputtering summaries indicate a sputtering ratio 
for each atom. Figure 123 shows the normalized sputtering 
summary for the randomly selected impact points. Figure 124 
shows the same information using the grid of impact voints. 
The four groups of 36 randomly selected impact points had a 
Mean Sputtering ratio of 0.812 with a variance of 0.0048. 
Additional runs ( not shown ), combined with the four 
Prevlous QGrOUpS to give a total of 21 groups of 36 randomly 
selected impact points ( 756 shots ), resulted in a mean 
Seer erang ratio of 0.782 with a variance of 0.016. fTfhis 
compares very favorably with a sputtering ratio of Q.778 
mong the grid Of impact points. Figures 125 and 126 show 
the "probability" results for the randomly selected impact 
points and the grid of impact points respectively. It is 
quite avident that they are almost identical. Figure 127 
illustrates the angular dispersion results for randomly 
Selected impact points. A definite shift of the tiignhter 
1sotopes toward the surface normal is observed. Figure 128 
illustrates the same information for the grid of impact 
Boints. The absence of the shift of the light isotopes 


toward the normal is noted. 
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moe OLMULATION MODEL 


The more sophisticated a simulation model becomes, the 
Mere Uncertain are its results. SOucmor cnese Nnceortaint los 


are now beginning to appear in this model. 


Pies LattIiCe “generators can now construct randonly 
positioned isotopes. Symmetrical unfolding of the event 
results is the reverse of the folding required to reduce the 
representative impact area. The model can now display the 
ejection angle of the sputtered isotopes in a histogram for 
Visual comparison. It now determines the distance from the 
surface layer center to each ring of sputtered atoms, the 
number of atoms and the number of times each atom sputters 
MmamedGch Cing, and the probability that that each ring 
contained an atom that was sputtered. Other modifications 
were made in the model to approximate more nearly 4 
meaningful physical decription. These include changes in 
the distances from the lattice surface at which calculations 
acre truncated, the depth of lattice needed to prevent any 
change in the yield, and a systematic determination of the 
distance above the surface at which an atom is declared to 


have sputtered. 
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Bee een lt aINMENT 


The containment of an event is paramount in a complete 
description of the sputtering process. However, it is less 


Significant when studying mechanisms and trends. 


Based on the definition of containment ( no increase in 
yield with an increase in the microcrystallite size ), the 
100-eV results showed that a 9x4x9, (100) orientation 
Microcrystallite contained only the helium and xenon events. 
The 100-eV, mercury event 1s contained in an 11xUxil 
Vote ce . To contain a 100-eV copper event, which had the 
largest yield, a 13x4x13 lattice was needed. Sputtering 
mechanism traces revealed the importance of boundary atoms 
in deflecting and channeling atoms in the Spuccer. 1g 
process. By choosing the PKA to be an atom other than the 
one at the surface layer center, the effect of an increased 


lattice size could be duplicated. 


Investigation with a 17x4x17 microcrystallite revealed 
that Gibson-II potential function simulations are more 
Guemeciit tO Contain than when using the KSE-B potential 
function. Only the 100-eV, copper event was contained using 
Gibson-II, while using KSE&-B the 100-eV and 200-eaV, copper 
events were contained. It is quite obvious that to contain 
any higher-energy events would require a microcrystallite so 


large that the computer run time would be prohibitive. 
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C. ANGULAR DISPERSION 


Because the low-energy events resulted in so fow 
Beuttered atoms, the angular dispersion results for any 
particular ion are inconclusive. A trend of light isotope 
enrichment is observed when the results from numerous runs 
are combined. For example, the histograms of the four sets 
of 36 randomly selected impact points show a definite shift 
Seecne lighter isotopes to snaller ejection angles. In his 
report on the enrichment of lighter isotopes in the normal 
direction, Wehner [24] noted only a 0.6% enrichment in Cu 
M3) . The number of sputtered isotopes was not sufficient 


to detect this small difference. 


Toe "reflective collision" theory, which Wehner suggests 
to be the cause, was observed to have very little effect on 
the sputtered atoms. Investigation revealed that in most of 
the simulation runs, the same atoms sputtered independent of 
its isotope type or the various positions of tne lattice 
isotopes. When the PKA was a light isotope, the atoms 
that did sputter were ejected more nearly normal to the 
Surface. Varying the isotope type of the PKA, atoms that 
were known to sputter, and some of the target atoms that 
they collided with resulted in differences in the ejection 


angle normal of less than 0.4°. 


De RING FORMATION 


Results dofinitely show that the sputtering mechanisas 
Momeprodice a formation of rings of sputtered atoms around 


the impact area. The results indicate that the transfer of 
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fementum along the <100> and <101> directions in the surface 
layer is relatively easy. The predominance of sputtering of 
next-nearest neighbors, in the <101> direction is evident. 
For example, at 100-eV, all simulations using an ion of the 
mass of argon or greater sputtered the same atoms, but with 
different frequencies. The particular atoms that sputter do 
not appear to be dependent on the mass of the incident ion. 
However, the number of times the atom sputters does appear 
to be dependent on the ion mass. Ring and atom sputtering 
probabilities clearly show the sputtered atoms from certain 
Site directions and the distances of the sites from the 
center of the PKA. 


A trend was noted: as the energy and mass of the ion is 
increased, the sputtered atoms come from sites which appear 
to spread outward from the impact area. These higher-energy 
events give the appearance of a "liguid splash", with the 


size of the crater dependent on the ion energy and mass. 


Pee Lon /7ATON MASS RATIO 


iemmesults of this investigation do not support a 
MOnatonic increase in the sputtering ratio with increased 
ion mass as implied in figures of Sigmund{[28] and Eer Nisse 
fee? ). All the plots showing the sputtering ratio as a 
Gunction of the ion mass, clearly indicate a maximum near an 
Mon/atom mass ratio of one. At ion/atom rnass ratios greater 
than one, the sputtering ratio appears to become asymtoptic 
momma. Lower value. This may be an artifact of the computer 
model. These plots also highlight the containment problen. 
Gmemeoehavior of the sputtering ratio, for larger ion/atonm 
/Missebatios, is erratic because the microcrystallite does 


not contain the event. 
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The general shape of the sputtering ratio curve, when 
plotted as a function of the ion/atom mass ratio, is in 
agreement with the experimental data of Wehner and his 
co-workers. Results conclusively show that the ion/aton 
mass ratio 1S an important factor in the determination of 


the sputtering ratio. 


PeeUANDOMILY SELECTED IMPACT POINTS 


The sputtering ratio obtained using the grid of impact 
points is extremely close to that obtained when randomly 
selected impact points were used.The similarity between the 
Sputtering summaries and the atom sputtering probabilities 
are quite evident. The quantity of data obtained with the 
randomly selected impact points was reguired to confira that 


Tight isotope enrichment was present. 


Based on the comparison of results obtained from the 
randomly selected impact points and from the grid of impact 
Pests, the systematic method of "walking" through the grid 


seems completely justified. 


Gee CULTURE RESEARCH 


Even the more than 1000 hours of computer tine used in 
running these simulations did not begin to encompass the 
permutations and combinations of parameters that need to be 


investigated. 
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The (110) and (111) orientations need to be studied in 


depth. The few sputtering mechanisms observed were? th2 same 
memeetose fOr the (100) oLlentation. However, the sputtering 
Mmat2o aS a function of the ion mass, containment 


requirements, ring formation, and angular dispersion of 
isotopes need further investigation. nie is highly 
recommended that the effect of momentum transfer along 


certain directions be investigated. 





APPENDIX A 


ISOTOPE GENERATION 


All atoms in the lattice were one of two types; type 1 ~ 
the heavier, least abundant isotope, or type 2 - the 
lighter, most abundant isotope. The average atomic weight 
of copper is 63.54 amu. The isotope weights are 62.93 anu 
and 64.93 amu with natural abundances of 69.1% and 30.9% 


respectively. 


The lattice was generated witn either the isotope masses 
and their corresponding abundance or the weighted average of 
the isotope masses with an abundance of 100%. Initially, 
all atoms were forced to the type 2 atom in the lattice 
generator. Tf the isotope lattice was desired, the 
International Mathematical and Statistical Library (ISML) 
routine GGUB was called. This routine generates a uniforaly 
distributed set of psuedo - random numbers between zero and 
one. Each of these random numbers waS asSigned toa 
particular atom and was tested to see if 1t was less than or 
equal to 0.309. If the number was less than this test value 
(the abundance of the least abundant tsotope), the atom was 
Changed to a type 14 atom. A seed for the random number 
generator, for the particular size lattice to be generatei, 
Was chosen by trial and error which would produce the 
correct number of atoms above and below the test value to 


assure the correct isotope ratio. 
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A feature was designed to allow any number of the atom 
Masses (types) to be selectively changed after the lattice 
had been generated. Migomeredeure | Was ~ancorperated to 
investigate the changes and effects of picking heavy/light, 
heavy/heavy, or light/light collision atoms. This would 
facilitate the dynamic investigation of the "reflective 


collision" theory proposed by Winters and Sigmund [10]. 


The dynamics of the simulation remained unchanged except 
that all calculations that were dependent oon target atom 


mass were modified to use the correct nass. 
The simulation was now capapdle of generating a lattice 


with the correct isotope ratio and masses with randomly 


Selected sites. 
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APPENDIX B 


LATTICE FOLDING 


As mentioned in the simulation model development, a 
lattice that was symmetrical about the folding axes was 
reguired. Once this was achieved by the lattice generators, 
it was necessary to fold the lattice ina manner such that, 
after the folds, the reduced representative impact ar2a was 
identical to the complete representative impact area. After 
these folds were performed correctly and in the right 
Sequence, the results were identical to that obtained 1f the 
complete representative lmpact area had been bombarded 


systematically. 


After each fold, the sputtering count for all atoms that 
were Symmetric about the folding axis were summed. The 
Summing process started with the first atom in a lay2r to 
the highest numbered atom in that layer After this process, 
the sputtering count for each atom was twice what the sum 
was Since the sputs were summed each time a symmetric atom 
was encountered. Dit sceeoun Eng umer fect mor Gach told was 
Borrected by scaling the sputtering count of each atom after 
all the folds in a layer were completed. All necessary 


folds were then performed on each successive layer. 


66 








A. fmeO) ORIENTATION FOLDING 


Figure Zz shows the (100) orientation reduc2d 
representative impact area folding geometry. The reduced 36 
impact point area, when folded, represents the complete 288 


impact point area. 


Eemce che (100) orientation impact area is a 459-4590 
triangle, three folds were required to complete the folding 
process. The first fold was along the diangon The first 
Mold was along the diagonal axis, the second along the 
horizontal axis, and the third along the vertical axis. The 
diagonal fold process for all layers was the same. However, 
the horizontal and vertical folds were modified depending on 


which layer was being folded. 


After all the folds in a layer for the (100) orientation 
were complete, the sputtering count for each atom was 


Geyaded by eight to correctly scale the results. 


B. (110) ORIENTATION FOLDING 


meaiceet2 also shows the (110) orientation reduced 
representative impact area folding geometry. The reduced 36 
impact point area, when folded, represents the complete 144 


impact point area. 


Since the (110) orientation reduced impact area is a 
rectangle, two folds were required to complete the folding 
process. The first fold was along the horizontal axis and 


the second fold was along the vertical axis. The diagonal 
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mermdewas NOt necessary. The horizontal and vertical folds 
(exactly the same as for the (100) orientation) were 


modified depending on which layer was being folded. 


Meect agit the folds in a layer for the (110) orientation 
were complete, the sputtering count for each atom was 


Mmevarded by four to correctly scale the results. 


ce. (111) ORIENTATION FOLDING 


Mitemegt!) Orientation was by far the most difficult to 
determine. Figure 13 shows the (111) orientation reduced 
representative impact area folding geometry. The reduced 21 
lmpact point area, when folded, represents the complete 126 


impact point area. 


Since the (111) orientation impact area is a 309-609 
triangle three folds along the diagonals resulted in the 
overlap shown in Figure 13. However, if the set of three 
folding sequences were summed the results could be scaled to 
represent the complete impact area. The number one fold was 
along the positive slope diagonal axis, the number two fold 
was along the negative slope diagonal axis, and the number 


three fold was along the vertical axis. 


Before each sequence of folds, the sputtering count for 
each atom was reinitialized to the original value. The 
first sequence (1,2,3) was performed and upon completion the 
Sputtering count for each atom was divided by eignt to 
correctly scale the results. The second sequence (1,2,1) 
was performed and again the sputtering count for each atom 
was divided by eight for scaling. The third sequence 
(1,3,2) was performed and again the sputtering count for 


each atom was divided by eight for scaling. 
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Finally, the results obtained from each sequence was 
summed and the final sputtering count for each atom was 
divided by four to correctly scale the results due to the 


Summing process. 
Since the geometry of the (111) Orientation was 


complicated and not fully developed, the folding process was 


performed only for the first layer. 
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GLOSSARY 


ABUNT, ABUN2 


The abundance of the target's isotopes. 


AC 
A distance measurement used in impact point 


generation. 


ALFA, ALFA2 
The cosine of the angle between vectors R1,R2, - 


squared. 


ALPHA 
Input Morse potential parameter. 
ANGLE 
The angle between the direction of motion of an 
ejected atom and the surface normal in degress. 
AVE 


The sputtering ratio. 


AXT, AYT, AZT 
Components of the random thermal displacement. 


BENRGY 


The lattice binding energy in ev. 
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BIND 


BMAS 


BQZ 


BUL 


BX, BZ 


BXMAX, 


BZMAX, 


CELS 


CGB1, 


CGD1, 


CGF1, 


Negative of the total potential energy (TPOT) at 
time Zero. Retained for energy conservation 


checks. 


The mass of the bullet in amu. 


MRCGmGeyYstal COnStant —({(1/ROOT(3)). 


The bullet heading Hollarith. 


woscalved xX,z coordinates of the impact pdoint. 


BXMIN 


X direction boundariss of the target impact area. 
BZMIN 


Z direction boundaries of the target imoact area. 


Menaareelonal force multiplier when frictional 


forces are included in the progran. 


CGB2 


Morse potential parameters. 


CGD2 


Morse potential parameters. 


CGF2 
Mopse torce function parameters. 
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imo, CFI, CF2 


Force parameters of the cubic fit between 


alamwsorn—-Mayer functions. 


CO, CORF 


"Phonon" displacement correlation parameter 


aS an option in WARM. 


CORN 
i CORT 


Bom, COY, COZ 


Morse 


used 


Directional cosines of the incoming ion's velocity 


vector, relative to the surface normal. 


Sort, COYI, cozi 
Negative values of COX, COY, COZ. 


20, CPi, CP2, CP3 


Potential parameters of cubic fit between Morse 
and Born-Mayer functions. 
CVBD 
A ratio used to avoid repeated division. 
CVD 
Distance conversion factor - angstroms to meters. 
CVE 
Energy conversion factor - ev to joules. 
CVED 


A ratio used to avoid repeated division, CVE/CVD, 


used to convert exponential potential functions to 


forces. 
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CVM 


Nemomconversion factor - amu to kg. 


CVR 


Conversions tactor —- lattice units to angstroms. 


DBXMAX, DBZMAX 


Upper boundary limits for the target impact area. 


DCON 
Input Morse potential parameter. 
DFF 
Time Gistance closer than ROE that an atom 1S to 
the bullet. Used when forces are turned on or off. 
BiLSTt 
The distance between any two atoms. 
DNN (I) 


Dstamee to the [lth ring in lattice units. 


DRX, DRY, DRZ 
Components of DIST. 


DT 
Length of each timestep in seconds. 

DTI 
Number of lattice units the most energetic atom 
Can move in one timestep. 

DTOD 


DT over distance - a ratio used to avoid repeated 


deeyi Sion. 
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DTOM1, DTOM2 


DTOMB 


Pipewer target 1Sotop]S mass - a ratio wsed to 


avoid repeated division. 


DT over bullet mass - a ratio used to avoid 


repeated division. 


DX(I), DY(I), DZ(I) 


EMAX 


Pres 


PEAR 


EV 


EVR 


EXA, EXB 


FA 


Suande = ==n the position of the Ith atom. 


The velocity sguared of the fastest atom. Used to 


make ETEST. 


The kinetic energy of the most energetic aton. 


Energy threshold to escape the surface. 


Primary energy of the bullet in ev. 


Primary energy of the bullet in kev. R indicates 


a read-in and saved value. 


Input Born-Mayer potential function for the target 


PeOMmmcalrgetatom interaction. 


The temporary component force increment on an 


acon. 
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PAC 


Pern SH 


FLAT 


Jj 


PMAX 


FORCE 


eed @ 


FRC 


FRAN 


FRANR 


tmemGccEement Size USed in SETUP in docating the 


Leonean pact point. 


The elapsed program running time in seconds. 


Miteerec Unit in angstroms. 


A small number used in checking potential energy 


ZeEO pont. 


The maximun frequency of occurrence in *he 


histogram intervals. 


Manet heat Value of the force Eunction for a 


SpcemeLte spain Of atoms. 


ERvemconnective £orce at ROS. 


Numerical value of the target force function at 
ROE. 


A real variable for number of passes conplezed at 


termination. 


The reciprocal of FRAN. 


Us 
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FREQ1, FREQ2 
Number of events in each inser va Or the 


histogran. 


FSCAL 
MIeStogran scale factor. 

FTDI 
A parameter used to determine DT by maximum energy 
method. 

FXA 


rd 
A Born-Mayer force function parameter. 


mote), FY(1), FZ(TI) 


Components of the total force on the Ith aton. 


HBMAS 
Halt Builet Mass - a ratio used to avold repeated 
dave Sion. 

HDTOD 
Hdttep bs Over “~dlLstance = da Tatio used to avoid 
repeated division. 

HDTOM 
Half DT over most abundant isotope mass - a ratio 
used to avoid repeated division. 

HDTOMB 


Half DT over bullet mass - a ratio used to avoid 


repeated division. 
HTMAS1, HTMAS2 


Half target isotope mass - a ratio used to avoid 


repeated division. 
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IBOX1, 


ICOUNT 


IDEEP 


ees 


1a ge 


IFLAG 


IFM 


cr T 


IH1 


I fel Me 


[BOX2 


Ejected atom angle tally keeper. 


Plag used in sequencing. 


The alc Sie frozen layer. aig the target 


microcrystal. 


Elapsed time in 2.6 microsecond blocks. Used to 


determine actual elapsed time. 


The atom number of any target atom that has 1ts 


mass constrained due its isotope ratio. 


Flag used in sequencing. 


The isotope type (heavy or light) of any target 
atom that has its mass constrained due to its 


lsotope ratio. 


A number used aS a constant when generating 


pseudo-random numbers. 


Alphanumeric array of the program heading. 


Alphanumeric array of the Morse parameters. 
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IHB(1), IHB(2), IHB(3), IHB(4), IHB(5), IHB(6) 


Alphanumeric array describing the bullet. 


IHBPOT 
Mhbeeeype of ion-atom potential. 
LTHBULL 
The bullet elemental symbol or name. 
IHS 
Alphanumeric array for the type and orientation of 
the crystal. 
in tT 


PeenemImeriG aceay describing the target. 


IHT(1), IHT (2), IHT(3), IHT (4), IHT(5), IHT(6) 


Alpha numeric array describing the target. 


IHTARG 
The target elemental symbol or name. 
LHTPOT 
imo motetakGetsacOle> target atom porential - 
ILATT 
A flag used in the lattice generators. 
ILAYER 
A variable used to identify which layer is being 
i Onded. 
ILL 


The number of the last atom in the first layer. 
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£MY 


IN 


INDEX 


‘INT 


EOR 


IOUT 


PEAS S 


IPNUM 


IPNOUMS 


IPROG 


IRAN 


The mass of each target atom with the correct 


isotope ratio. 


The number of intervals in the histograms. 


A variable used in determining phase in a dynamics 


cycle step. 


The width of the intervals in the histograms. 


ate 


The target plane of incidence. 


Flag used in sequencing. 


The number of passes (shots) completed. 


The impact point number. 


Scaled impact point number. 


Input variable used to skip attractive potential. 


Random number seed used in NORMAL subroutine which 


generates pseudo-random numbers. 


US, 





IRANGE 


ivemuppenr Limlt in the histograms. 


IREAD 
The number of atoms which have the isotop2 type 


constrained externally. 


ISEEDA, ISEEDB, ISEEDC 
Seeds used in the uniform random number generator 
which determines the isotope ratio and 


Gest rnibuetLOnN . 


FSET 
Flag used in sequencing. 
ISHUT 
A parameter used to shut down the progran. 
ISNUM 
The atom number of a sputtered aton. 
ISOGEN 
A flag used as an option to generat2 the target 
atcm masses with 1sotope or single element values. 
ISPX, ISPZ 
Pio tiaiestdrt POlNt x,Z coordinates. 
Jhsk 
Unscaled x coordinate used in lattice generation. 
Jt! oO 
The lattice temperature in degrees Kelvin. 
LEELA E 


Flag used in sequencing. 
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Itt 


An odd-even integer used to determine atom site 


Location. 


Ix, LY, 12 


Pie saumoer OL X,y,z, Planes in the crystal, 


mar, LYP, IZP 


Seyccal Gimensions in X,Y, Z- 


JN 
Biestimertcalvvyalue Of tne histogram int 2rva ls, 
JOUT 
A histogram sequencing variable. 
JSCAL 
A histogram scale factor. 
JT 
Unscaled y coordinate used in lattice generation. 
Bp) 
Variable used to establish atom site. 
weT 
Variable used to establish atom Site. 
K2SUM 
The square of the sputtering count which is used 
in calculating the variances. 
KATOM 


The number of atoms in a ring. 
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ReeRST 


KRING 


KSPUT 


KSUM 


KT 


LCUT (I) 


LD 


LL 


HeZ, LL3 


LNOUM (I) 


LNUMF (1) 


A sequencing variable. 


The highest numbered ring in a layer. 


PiomiIioecr OL Sputs £rom all atoms in a ring. 


The total sputtering count before folding. 


Unscaled z coordinate used in lattice generation. 


Used to identify an Ith atom which can be omittad 
EE0m the force and potential calculations oby 


Mime eLeULr( i) fron Zero to one. 


The highest numbered atom in a mobile layer. 


The highest numbered atom in the entire crystal. 


turecesand triple the value of LL. 


The number of times the Ith atom sputters aiter 


Sm esnOus in the target impact area- 


The number of times the Ith atom Sputters after 


all shots and the summary is complete and folded. 
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LNUMI (I) 


JES) 


bas) 


LSUR 


LTYPE 


MAX 


MCRO 


NY 


NCT 


NCTT, ND 


The number of times the Ith atom sputters after 


all shots but before folding the summary. 


Sum of the Miller index integers. 


Visine used to identify the type of surtace. 


The lattice surface type. 


Ticmiyoe Or lattice {1.e. fec ). 


A height scale factor used in the histograms. 


The number of simultaneous eguations solved by 
CROs. 4. 


NLAY - 1, a sorting variable. 


Increments the timestep interval. 
a 


Data output interval in number of timesteps. Also 
the timestep interval between potential energy 


Galeulatvons. 


Q2 








NF 


A sequencing variable. 


NFIRST 
| The number of the first atom in a given layer. 
NLAY 
The highest numbered atom in a given layer. 
NLAY1 
The highest numbered atom in the first layer. 
NLINE 
Rice line number in the outout block. 
NOTH 
A Hollarith consisting of blanks. 
NPAGE 
A page numbering variable. 
NRAN 
Variable for number of passes (shots) completed at 
termination. 
NRS 
A restart flag. 
NS 


The next timestep number at which output will be 


Dern ted . 


BY 





NSHOTS 


The number of passes (shots) to be completed in a 
program. Used to break up an extremely large, 
many-shot program into a few small programs with 


manageable computer run time. 


NSPUT (T) 
The number of atoms that sputter for tne Ith 
impact point. 
NSR 
The first timestep where output 1s printed. 
NSS 
Pm derompraie FLfal Output status of curpent shot. 
NT 
The timestep number. 
NT 
Timestep number limit before shut down. 
PAC 
Parameter for bullet force function correction. 
PBMAS 
Bullet mass in kg. 
PEXA, PEXB 
Input Born-Mayer potential function parameters for 
the bullet-target interaction. 
PFIV 


Bomndt  tasye. (0.5). 
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PFPICc 


PFRC 


PFXA 


PKE (I) 


PNEG 


PNUM 


POT 


PPE(I) 


BeOS 


PPTEC 


PROB 


puleMet wEOrce function evaluated at ROE. 


Bullet-aton 
at ROE. 


repulsive force interaction evaluated 


Bullet force function parameter. 


The kinetic energy of the Ith atom. 


Total negative potential energy in the lattice. 


Alphanumeric variable for the impact point number. 


The potential energy between two atoms. 


The potential energy of the Ith aton. 


Total positive potential energy in the lattice. 


Bullet potential function evaluated at ROS. 


treme probability that an aton will sputter from a 


given ring. 
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BC 


PTE (fT) 


QM 


QUIT 


R1, RIS 


R1, 


R2 


R3 


Ra, RZ 


R2SQ 


R1Y, 


Target atom potential function evaluated at ROE. 


The total energy of the Ith aton. 


A small number used in checking the kinetic é4energy 


Zeno POl nt. 


Cocenre variable checked against the total 
petential energy which can be used as a shut down 
eEiterion. 

Diseamce Etrom (0,0) impact point to the initial 


Bullet position. 


R1Z 


Pie x,y, 2 COMNpOnents of RIS. 


Magnitude of vector from impact point to the first 
meonenrt by the bullet. 


Maommtude Of vector from impact point to buliest 


start position. 


SGalar product of vectors RT, R2. 


R2 squared. 
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RAD 
The angle between the direction of an ejected atom 


and the surface normal in radians. 


RADDEG 
Conversion factor - radians to degrees. 
RAN (I) 
A list of normally distributed random numbers. 
RANGE 
A scaling variable used in the histogran. 
RBX, RBZ 
Unscaled x,z coordinates of the impact area 
reference point. 
RDLST 
The distance between an atom and the center point 
of its layer. 
RE 
Input Morse potential parameter. 
RINCR 
The RANGE increments. 
RMIDPT 
The midpoint of each interval in the histogran. 
BOE. nOEZ 
Nearest neighbor distance - squared. 
ROEA 


The maximum cutoff for the Born-Mayer potential. 
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ROEB 


The minimum cutoff for the Morse potential. 


ROEC, ROEC2 
The maximum cutoff for the Morse potential - 
squared. 

ROEM 
Defines the region where the repulsive force 
function must be modified. 

mexett), RY(I), R2Z(I) 
Wye caecooLrdinates of the Eth atom at any time. 

RXBL, RYBL, RZBL 
Maximum negative x,y,Z2 coordinates an atom nay 
Move before it is LCOT. 

Rese, RYBA, RZBM 
Maximum positive x,y,z coordinates an atom may 
move before it is LCUT. 

ae, Ril, RZF 
Final x,y,z coordinates of the Ith atom when it is 
eur. 

MEE(I), RYI(1), RZI(TL) 
tiemieneaton's initial x,y,z coordinates. 

BAR (LE), RYK (7), RZK (T) 
Temporary X,y,Z coordinates of the Ith atom during 
the force cycle. 

Ba&ay  kYS, RLS 


Scaled x,y,z coordinates of the impact point. 


39 





SAVE 


SBX, SBZ 


One-half POT. 


Scaled x,z coordinates of the impact point. 


SCALE 
A scaling variable for the frequency of events in 
the histogran. 

pew, SCY, SCZ 
The x,y,z lattice scale factors used to convert 
the generated lattice to match a real crystal. 

sexR, SCZR 
The inverse of SCX and SCZ - a ratio used to avoid 
repeated division. 

Sp LST 
The sguare of the distance from the center of a 
toaven tOudn AtOMm. 

Seep Be 
A library subroutine used to keep track of run 
eames 

SLOW 
A cutoff velocity variable checked against a long 
Bl. 

SPDX, SPDZ 


Meme cEaction of the target impact area that the 


impact points are incremented. 
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SPX, SPZ 


The x,Z Components of the distance from the impact 


Seca tererence point to the impact point. 


SSCZ 
feeeeeeesoeake factor used in the (111) “orientation 
lattice generation. 
SLD 
The standard deviation of the sputtering ratio. 
 STOPER 
This variable stops calculating and saves 60 
Seconds fOr COMmpillation and output. 
SUR 
Alphanumeric variable describing the surface. 
TA, TAMP 
The thermal amplitude of the lattice. 
TAR 
Alphanumeric variable describing the target 
element. 
TARGET 
Alphanumeric variable describing the target 
element. 
TAS 
A scale factor - TAMP/ROOT (3.0). 
TASF, TASN 


Random thermal amplitude scale constants. 


on 





TE 


TEMP 


west 


ar AC 


THERM 


TIME 


ini St, 


meASl, 


TPKE 


2PoT 


TRUN 


Pilemccral Kinetic and potential energy of all the 


atoms. 


The temperature of the lattice in degrees Kelvin. 


- 


A termination test of energy. 


A time factor ratio used to determine DT by 


maximum force methods. 


Mean thermal energy of a lattice atom. Used in 


Shut-down tests. 


The elapsed calculating time in seconds. 


Tile Z 
The titles for the histogram outputs. 


TMAS2 


The target isotope masses in amu. 


Tremtotai Kinetic energy of all the atoms. 


The total potential energy of all the atoms. 


The time the program is to run in minutes. 


gZ 





VEL 
DHemvetocity calculation - SORT (ev/HBMAS) . 


¥Sss 
Meecetiporary Storage Variable for the velocity 
components. 

me EST 
A test velocity used to determine if the atom has 
sufficient velocity (energy) to escape the 
surface. 

WIS 


A scale factor - 0.8*VTEST. 


Wet), VY(I), VZ(I) 


Meek, ¥,Z COMpONnents of the Ith atom's velocity. 


VXS, VYS, VZS 
Gime ye2 ee COMDONeEeNtS of =the bullet's tnitial 


velocity. 


ZDIST 
The x component of the distance from the center of 
ameaviels EO the atom. 

XMAX, XMIN 
ivesx planes to determine whether two atoms are 
close enough to interact. 

YLAX (T) 


Piomersystal relaxation of the [th layer in the ~y 


a 


Geercece 1oOnNn. 
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YMAX, YMIN 
The y planes to determine whether two atoms are 


close enough to interact. 


ZDIST 
The z component of the distance from th2 center of 
a layer to the aton. 

ZE 
Absolute zero (0.0). 

ZMAX, ZMIN 
The z planes to determine whether two atoms are 
close enough to interact. 

moO P 


Used to print out a warning header that a pass was 


terminated before completion. 
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Figure 54 - Angular Dispersion - (190), 1-k2aV¥, Cu/Cut 
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